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Making use of the Bogoliubov-de Gennes equation for the d-wave superconductors, we investigate 
the quasi-particle spectrum around a single vortex. Taking pp£ = 10, we found that there are bound 
states which are localized around the vortex core, and extended states which are rather uniform, for 
\E\ < A where E is the quasi-particle energy and A is the asymptotic value of the order parameter 
for away from the vortex. 



Much attention has been paid to the quasi-particle 
spectrum around a single vortex in d-wave superconduc- 
tors, since the discovery of the high-T c cuprate supercon- 
ductors Q-^j. As is well known, d-wave superconductiv- 
ity in both the hole-doped and the electron doped high-T c 
cuprates has been established ||-^]. 

Recent scanning tunneling spectroscopy (STS) sees dif- 
ferent core structures of the vortex in YBCO and Bi2212 
|j] |llj. Contrary to the theoretical expectation a) 
they observed a bound state with energy E, a fraction of 
A, and b) they have not seen any clear fourfold symme- 
try. 

Motivated by these experiments, Morita et al. (T^| have 
abondoned the quasi-classical approximation used in ear- 
lier analysis B and proposed to study the bound states 



in terms of Bogoliubov-de Gennes equation |13|. Indeed, 
this approach appeared to give the correct description 
of the observation. However, later Franz and Tesanovic 
claimed that there should be no bound states HK Fur- 
ther this claim was confirmed by Yasui and Kita |l5| and 
by Takigawa et al. | fitfl later. In a vortex in s-wave su- 
perconductors, Caroli et al. JH| have shown there are a 
series of bound states. Further a detailed structure of the 
bound state wave function is explored later by Gygi and 
Schliiter [fL8[ . As we have shown later, there are many 
bound states in d-wave superconductors as in s-wave su- 
perconductors We don't know for sure the origin of 
this disagreement. But the possible source of their er- 
rors in previous studies |l^-|l6| is easy to locate. They 
neglected in their analysis the conservation of the angu- 
lar momentum around the vortex. Of course the strict 
conservation of the angular momentum is broken due to 
the fourfold symmetry of A (fc). On the other hand the 
angular momentum is still conserved by modulo 4, and 
this is adequate to gurantee the presence of bound states. 
Therefore the structure of the bound states in the vicin- 
ity of the vortex core of d-wave superconductor appear 
to be very similar to the one in s-wave superconductor. 



Then the only clear difference is the presence of extended 
states first discovered in |l2j . Also these extended states 
give rise to the Volovik effect [|l9|, which has been ob- 
served as the V^ff-term in the specific heat J2(| and more 
recently as the H linear term in the thermal conductivity 

Earlier we have considered extreme quantum limit, i.e. 
Pf^ ~ 1 [|l7]], but recent experiment shows pfS, — 20, 30 
for Bi2212 and YBCO, respectively Therefore, in 

this letter, we take pp£, = 10. We consider quasi-two 
dimensional d x 2_ y 2-w&ve superconductors in a disk with 
radius R. We put a single vortex at the center of the 
disk. When R is large enough, the magnetic field is weak 
and the vector potential is small, so we ignore the vec- 
tor potential. The Bogoliubov-de Gennes equations are 
given as, 
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d y A (r) d y ] u n (r) = E n v n (r) . (2) 
The order parameter is given as, 

A(r)=Y,^T(l-V(En)) 

n P F 

x [(d x u n (r)) (d x v* (r)) - (d y v n (r)) (d y u* n (r))] . (3) 

In the following we assume that the vector potential takes 
following form, 



A(r) = |A(r)|e 
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We use the Fourier-Bessel expansion following Gygi- 
Schliiter as, 
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where </> TO j is Fourier-Bessel basis given as, 



</>mj 0) = 
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Here a m j is the j-th zero of the Bessel function of the 
m-th order. Then the Bogoliubov-de Gennes equation 
becomes as, 
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A^ m _ 1 's are matrix elements and given as, 
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where <j> d mj (r) = RJm ^ (oim]) J m -i and 0*^. (r) = 

flj — -jJm+i ( a "i{ r )- The order parameter is now 
given by, 

n jij 2 rn 

[ftmji 0) 0t+3 3 2 ( r ) Unm 3l V nm +3 32 
+4>mh ( r ) <t>m-\n 0) UnmhVnm-lh] ■ (12) 

In order to fix the particle density, we determine the 
chemical potential by the condition, 
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where N e is total particle number. 

From Eqs. ^ and ^|, there are coupled sequences of u m 
and v m and there are four sectors of them. 

0th • • • , ti n _4j , V n -ij , U n 0j , W„3j , • • • 
1st ' ' • , V n — 4j , U n _3j , WnOj , 1*nl j , ' ' ' 

2nd • • • ,v n -3j,u n - 2 j,v n ij,u„2j, ■ ■ ■ 

3rd • ■ • (14) 



There is a following symmetry between these sectors, 



Therefore we solve Eqs. (||) and (^|) for the 1st and the 
2nd sectors and using the symmetry we determine A (r) . 

For the numerical calculation, we take R = 15£, 
A Q /E C = 0.4 and p F t = 10. 

First, we show the temperature dependence of the or- 
der parameter in Fig. |TJ. 
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FIG. 1. Temperature dependence of the order parame- 
ter, (a) T/T c = 0.1, (b) T/T c = 0.2, (c) T/T c = 0.3, (d) 
T/T c = 0.5, (e) T/T c = 0.7, and (f) T/T c = 0.9. 



There is an oscillation of the order parameter away 
from the vortex core. This comes from the geomet- 
rical resonance of the system and the boundary con- 
dition. Comparing with the s-wave superconductors 
|l8|,P3|, there is no structure inside of the core. 
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FIG. 2. Local density of states as a function of E and r at 
T/Tc — 0.1. It is normalized by density of states of normal 
state N . 
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FIG. 3. Weight of bound states and extended states for 
E/A — 0.05 ± 0.02. 
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FIG. 4. Wave functions of a bound state, (a) |u(r)| and 
(b) \v (r)\ 2 . Its energy eigenvalue is E/A = 0.0481. 

We show the local density of states for T/T c = 0.1 in 
Fig. |. There is a small angular dependence, but we in- 
tegrate out the angular dependence for clarity. Because 
of the nodes of the d x 2_ y 2 order parameter, there is a 
finite density of states inside of the energy gap. Also 
there is a peak around r — and E = and a broad 
peak close to the boundary. Our previous calculation 
shows that for ppt; = 1.0, the peak around a vortex core 
moves to gap edge [ [l7| . In order to examine the origin 
of this peak, we show the relative weight of the bound 
state versus the extended state for E/A — 0.05 ± 0.02 
as function of r/£ in Fig. |^. The relative weight means 

a sum of J Q 27r u (r) 2 + v (r) 2 d6. As is readily seen the 



bound states are localized near the vortex core, while 
the extended states are rather flat. So this should give a 
quasi-uniform density of states at E = ||. These are 
shown is Fig. [| and Fig. |[ respectively. 




FIG. 5. Wave functions of an extended state, (a) |w(r)| 2 
and (b) \v{r)\ 2 . Its energy eigenvalue is E/A = 0.0331. 

Strictly speaking, we find other set of solutions, which 
cannot be characterized as the bound state or the ex- 
tended state. But we believe that they are due to the 
finite size of our disk which we have studied. 

The angular dependence of the density of states is 
shown in Fig. ||. There is an oscillation of the LDOS 
along the diagonal direction (x = y). This oscillation 
becomes large for large E. 

In Summary, making use of the Bogoliubov-de Gennes 
equation for d-wave superconductors JH|, we have stud- 
ied the quasi-particle spectrum around a single vortex 
with a weak-coupling model and for ppt; — 10. First, 
we find many bound states centered around the vortex 
core as shown in Fig. 0. This picture is very similar to 
the one found for s-wave superconductors i jTsfl . Second, 
ther are extended states with almost uniform amplitude 
near the vortex core (see Fig. |]). Third, there appear 
to be a small group of the mixture of the bound state 
and the extended state. However, we believe this mixing 
is the finite size effect. In order to clarify this question, 
we need a pararell study of the quasi-particle spectrum 
in disks with different size. There remains one mystery. 
Why STM in YBCO @ and Bi2212 |l) picked up only 
a single bound state? Perhaps some of bound states are 



3 



more readily accessible to STM? 




FIG. 6. Angular dependence of the local density of states, 
(a) E = 0, (b) E/A = 0.5, and (c) E/A = 0.7. They are 
normalized by the density of states of normal state No. 
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